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ABSTRACT 

This  investigation  has  been  tmdertaken  as  part  of  a  program  of  work  having  the  aim  of 
determining  a  suitable  fatigue  life  enhancement  option  for  the  non-circular  fuel  flow 
vent  hole  #  13  in  the  wing  pivot  fitting  of  the  F-lllC  aircraft.  The  stress  analysis  has 
been  imdertaken  using  two  dimensional  finite  element  methods  assuming  either 
elastic-perfectly  plastic  or  strain  hardening  material  behaviour  models.  The  effect  on 
critical  plate  stresses  due  to  enhancement  by  cold  expansion  and  interference  fitting, 
used  separately  or  together,  has  been  quantified  in  the  presence  of  representative  cold 
proof  test  loading  and  a  sample  remote  spectnim  loading.  For  aU  enhancement  cases 
considered,  irrespective  of  material  model,  the  stress  response  during  the  sample 
remote  loading  was  linear,  with  a  cycHc  stress  concentration  factor  of  unity. 
Enhancement  through  combined  cold  expansion  and  interference  fitting  was 
considered  a  great  deal  better  than  interference  fitting  only.  For  example,  the  combined 
enhancement  of  1.5%  cold  expansion  followed  by  0.5%  interference  fitting,  as 
compared  to  0.5%  interference  fitting  only,  led  to  a  change  in  critical  hoop  stress  from 
350  MPa  to  -1080  MPa.  These  favourable  results  indicate  that  such  an  enhancement 
procedure  would  potentially  be  suitable  for  extending  the  fatigue  life  of  the  F-111  #13 
fuel  flow  vent  hole  region,  pending  the  results  of  an  investigation  on  a  non-circular 
hole  geometry  currently  in  progress. 
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Stress  Analysis  of  a  Plate  Containing  a  Round 
Hole  with  Combined  Cold  Expansion  and 
Interference  Fitting  Under  F-lllC 
Representative  Loading  Conditions 


Executive  Summary 

An  area  of  major  concern  for  the  F-111  airframe  in  service  with  the  RAAF  is  the  wing 
pivot  fitting  in  which  there  are  a  number  of  non-circular  (elongated)  machined  fuel 
flow  vent  holes.  There  have  been  numerous  incidents  of  fatigue  cracks  at  fuel  flow 
vent  hole  #13  in  the  wing  pivot  fitting,  and  the  problem  could  compromise  the 
structural  integrity  of  the  F-111  fleet  out  to  the  planned  withdrawal  date  of 
approximately  2020.  Currently,  this  problem  is  being  managed  by  reshaping 
periodically,  the  fuel  flow  vent  holes  to  one  of  a  family  of  progressively  larger  shapes, 
thereby  removing  any  existing  cracks.  Unfortunately,  this  process  does  not  completely 
eliminate  further  cracking,  and  at  the  current  rate  of  crack  growth  and  hole  reshaping, 
this  method  may  not  be  sufficient  to  enable  these  aircraft  to  reach  their  desired  service 
life. 

This  present  investigation  has  been  rmdertaken  as  part  of  a  program  of  work  having 
the  aim  of  determining  a  suitable  fatigue  life  enhancement  option  for  the  non-circular 
fuel  flow  vent  hole  #  13  in  the  wing  pivot  fitting  of  the  F-lllC  aircraft.  The  stress 
analysis  has  been  undertaken  using  two  dimensional  finite  element  methods  assuming 
either  elastic-perfectly  plastic  or  strain  hardening  material  behaviour  models.  The 
effect  on  critical  plate  stresses  due  to  enhancement  by  cold  expansion  and  interference 
fitting,  used  separately  or  together,  has  been  quantified  in  the  presence  of 
representative  cold  proof  test  loading  and  a  sample  remote  spectrum  loading.  For  all 
eiihancement  cases  considered,  irrespective  of  material  model,  the  stress  response 
during  the  sample  remote  loading  was  linear,  with  a  cyclic  stress  concentration  factor 
of  unity.  Enhancement  through  combined  cold  expansion  and  interference  fitting  was 
considered  a  great  deal  better  than  interference  fitting  only.  For  example,  the  combined 
enhancement  of  1.5%  cold  expansion  followed  by  0.5%  interference  fitting,  as 
compared  to  0.5%  interference  fitting  only,  led  to  a  change  in  critical  hoop  stress  from 
350  MPa  to  -1080  MPa. 

These  favourable  results  indicate  that  such  an  enhancement  procedure  would 
potentially  be  suitable  for  extending  the  fatigue  life  of  the  F-111  #13  fuel  flow  vent  hole 
region,  pending  the  results  of  an  investigation  on  a  non-circular  hole  geometry 
currently  in  progress. 
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1.  Introduction 


An  area  of  major  concern  in  the  F-lllC  airframe  is  the  wing  pivot  fitting  manufactured 
from  D6ac  steel,  in  which  there  are  a  number  of  machined  fuel  flow  vent  holes  (FFVH), 
see  Figure  1.  Under  cold  proof  load  tests,  the  material  around  one  of  these  holes, 
FFVH13,  experiences  extensive  plastic  deformation,  which  results  in  the  introduction 
of  tensile  residual  stresses.  These  residual  stresses  coupled  with  the  local  material 
response  due  to  the  remote  loading  sequence  are  detrimental,  and  contribute  to  crack 
initiation  and  growth  at  the  affected  area.  There  have  been  numerous  incidents  of 
fatigue  cracking  at  FFVH13  in  the  wing  pivot  fitting  in  the  RAAF's  F-lllC  fleet,  and 
the  problem  could  compromise  the  structural  integrity  of  these  aircraft  out  to  the 
planned  withdrawal  date  of  2020. 


Figure  1:  F-111  Aircraft  &  wing,  showing  location  of  critical  hole  (FFVFI13) 


Currently,  the  problem  is  being  managed  by  reworking  the  elongated  fuel  flow  vent 
holes  to  a  family  of  progressively  larger  shapes  [1].  This  process  removes  small  cracks 
and  corrosion  as  detected.  The  extent  of  the  rework  depends  on  the  size  of  the  detected 
crack.  Unfortunately,  the  reworking  does  not  completely  eliminate  fxuther  cracking, 
and  at  the  current  rate  of  crack  growth  and  reworking,  this  method  may  not  be 
sufficient  to  enable  the  aircraft  to  reach  its  desired  service  life.  A  durability  and 
damage  tolerance  analysis  (DADTA)  is  presently  being  used  to  certify  die  use  of  these 
rework  shapes  as  a  valid  means  for  managing  this  critical  location.  If  the  reworking 
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procedure  cannot  allow  the  aircraft  to  reach  the  desired  service  life,  an  alternative 
method  of  life  extension  will  be  required.  In  view  of  this  problem,  AMRL  has  been 
tasked  with  the  development  of  a  proposed  non-circular  cold  expansion/ interference 
fit  plug  option,  for  the  life  extension  of  FFVH13,  with  the  aim  of  eliminating  crack 
growth  or  significantly  reducing  the  crack  growth  rate. 


The  typical  benefit  of  cold  expansion  and/or  interference  fitting  for  the  life  extension 
of  plates  containing  circular  holes  is  reasonably  well  established  and  documented.  For 
example,  the  use  of  interference  fitting  and/or  cold  expansion  has  been  reviewed  by 
Mann  and  Jost  [2].  Closed  form  theoretical  solutions  only  apply  to  the  simplest 
configmations,  so  typically  finite  element  analyses  are  employed  to  solve  problems  of 
practical  importance.  However  previous  work  does  not  give  clear  guidelines  for  either: 
(i)  non-circTilar  geometries,  or  (ii)  complex  remote  loading  conditions  and  realistic 
strain  hardening  material  properties.  Both  of  these  issues  need  to  be  addressed  when 
developing  a  life  extension  option  for  the  FFVH13  on  the  F-111  aircraft. 


In  this  report  we  address  the  second  issue  for  a  circular  hole  case,  as  a  precursor  to  the 
analysis  of  the  more  complex  non-circular  case.  The  key  parameters  in  this 
investigation  have  been  chosen  to  be  representative  of  the  Wing  Pivot  Fitting  region  of 
the  F-lllC  aircraft.  In  Section  2  the  process  definitions  for  cold  expansion  and 
interference  fitting  are  given,  and  their  effect  on  stresses  at  the  hole  edge  is  explained 
genericaUy.  The  faiite  element  methodology  employed  in  this  report  is  then  presented 
in  Section  3.  The  results  of  benchmark  analyses  are  then  given  in  Section  4,  where  a 
comparison  between  closed-form  theoretical  and  elastic-perfectly  plastic  finite  element 
analyses  is  made,  for  annuli  having  material  properties  representative  of  aluminium 
alloy  or  steel.  In  Section  5  the  problem  of  a  rectangular  plate  containing  a  circular  hole 
that  is  cold  expanded  and  interference  fitted  is  examined.  Here  both  elastic-perfectly 
plastic  and  strain  hardening  material  property  cases  are  considered,  and  the  plate 
stress  response  xmder  remote  CPLT  and  sample  remote  spectnun  loading  is  presented. 


2.  Cold  Expansion  and  Interference  Fitting  Process 

Definitions 


For  the  case  of  a  plate  containing  an  enhanced  circular  hole  it  is  helpful  to  highlight  the 
major  differences  between  cold  expansion  and  interference-fitting,  (see  also  Broek  [3] 
and  Jost  [4,5])  as  follows. 

In  the  Cold  expansion  (CE)  process,  die  hole  is  initially  expanded  to  a  level  sufficient  to 
cause  local  yielding  to  occur.  This  is  typically  achieved  by  passage  of  an  oversize 
tapered  mandrel  or  mandrel/  sleeve  combination  through  the  hole.  Then  the  mandrel 
is  removed  from  the  hole,  and  the  surrounding  elastically  deformed  material  forces  a 
reduction  in  hole  diameter  from  the  fully  expanded  size.  This  results  in  a  zone  around 
the  hole  containing  residual  compressive  hoop  stresses.  As  the  extent  of  yielding 


2 


DSTO-TR-0523 


during  expansion  is  increased,  the  larger  is  the  zone  of  induced  compressive  stress. 
Depending  on  the  degree  of  expansion,  compressive  reyielding  on  mandrel  removal 
can  also  occur,  although  the  size  of  this  reyielded  zone  is  substantially  less  than  that 
generated  by  the  mandrel  enlarging  the  hole.  Under  the  influence  of  cyclic  remote 
loading,  the  presence  of  the  compressive  zone  at  the  hole  boundary  results  in  the  mean 
value  of  the  induced  local  cyclic  stresses  to  be  significantly  less  than  for  a  non-cold 
expanded  hole.  Although  the  magnitude  of  the  local  cyclic  stress  range  is  typically 
unchanged  as  compared  to  the  non  cold  expanded  hole,  the  tensile  part  of  the  cycle  is 
substantially  reduced.  Hence  the  reduction  in  mean  local  stress  at  the  hole  botmdary  is 
typically  highly  beneficial  in  delaying  the  onset  of  crack  initiation  and  therefore 
increasing  the  fatigue  life  of  a  component  with  a  cold  expanded  hole. 

Interference  fitting  (IF)  is  the  process  of  installing  an  oversized  mandrel  into  the  hole, 
often  achieved  by  passage  of  a  tapered  mandrel  into  the  hole  until  a  required 
interference  level  is  obtained.  Unlike  the  cold  expansion  process,  the  mandrel  remains 
in  situ  and  typically  plate  deformation  around  the  hole  due  to  the  interference  fitting 
process  is  elastic.  This  process  has  two  effects.  Firstly,  a  residual  tensile  circumferential 
stress  is  induced  in  the  plate  at  the  hole  botmdary.  Secondly,  the  in  situ  mandrel 
provides  an  alternative  load  path  when  the  plate  is  subjected  to  remote  cyclic  loading, 
thereby  reducing  the  magnitude  of  local  cyclic  stresses  in  the  plate  at  the  hole 
botmdary.  Typically,  the  reduction  in  the  magnitude  of  the  local  cyclic  stresses  is 
highly  beneficial  in  delaying  the  onset  of  crack  initiation  and  hence  increasing  fatigue 
life.  However  caution  must  be  exercised  when  assessing  the  potential  benefit  of  this 
approach,  particularly  if  the  remote  loading  sequence  is  dominated  by  compressive 
loads.  In  such  a  case,  even  though  the  magnitude  of  local  cyclic  stress  at  the  hole  edge 
is  reduced,  the  higher  mean  stress  due  interference  fitting  signifies  that  both  applied 
tensile  loading  and  compressive  loading  will  be  fatigue  damaging.  Hence  the  potential 
improvement  in  fatigue  life  of  a  component  with  an  interference  fitted  hole  will 
depending  on  the  particular  remote  loading  sequence,  and  on  the  magnitude  of  the 
mean  stress  due  to  interference  fitting. 

Based  on  the  above  preceding  descriptions,  it  is  expected  that  greatest  benefit  to 
fatigue  life  extension  wiU  result  from  a  combination  of  cold  expansion  and  interference 
fitting. 


3.  Finite  Element  Method 


In  this  work,  elastic  and  plastic  plane  strain  analyses  have  been  xmdertaken  using  the 
PAFEC  finite  element  code,  level  8  running  on  a  Hewlett  Packard  K  series  9000 
computer  at  AMRL.  Eight-noded  isoparametric  quadrilateral  elements  (PAFEC  type 
36210)  were  used  for  aU  analyses  The  plasticity  subroutines  in  PAFEC  make  use  of  the 
Prandtl-Reuss  equations  in  conjtmction  with  the  von  Mises  yield  criterion  and  an 
isotropic  hardening  model.  For  all  plastic  analyses  loads  were  applied  incrementally, 
with  further  iterations  imdertaken  for  each  increment  until  acceptable  convergence  to 
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the  material  non-linear  stress-strain  constitutive  response  curve  was  achieved.  The 
convergence  criteria  was  judged  to  be  achieved  when  the  internal  energy  was  less  than 
1%  different  to  the  external  work  due  to  the  applied  loads,  for  a  given  increment. 


3.1  Geometries  and  Typical  Meshes 

The  notation  and  geometry  relevant  to  the  two  dimensional  idealisation  of  the  cold 
expansion  of  an  annulus  is  given  in  Figure  2.  The  outer  radius  of  the  annulus  is 
denoted  as  b  and  its  initial  iimer  radius  is  denoted  as  a.  An  oversize  mandrel  of  radius 
rm  is  inserted  in  the  hole.  The  resultant  hole  expansion  causes  local  yielding  in  the 
annulus  to  occur,  and  the  extent  of  yielding  is  defined  by  the  elastic  plastic  boimdary  c. 
Upon  removal  of  the  mandrel  from  the  hole  the  surrotmding  elastically  deformed 
material  forces  a  reduction  in  hole  diameter  from  the  fully  expanded  size.  Depending 
on  the  degree  of  expansion,  compressive  reyielding  on  mandrel  removal  can  also 
occur,  and  the  size  of  this  reyielded  zone  (which  is  substantially  less  than  that 
generated  by  the  mandrel  enlarging  the  hole)  is  denoted  as  p. 


Figure  2:  Notation  and  geometry  relevant  to  annulus  undergoing  cold  expansion. 
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Figure  3  shows  tihe  mesh  used  for  the  finite  element  analysis  of  the  annulus,  which  is 
not  subjected  to  remote  loading.  Due  to  the  axi-symmetric  nature  of  the  annulus,  and 
with  the  apphcation  of  appropriate  bovmdary  conditions,  it  was  only  necessary  to 
model  a  small  (2.5°)  segment  of  the  original  geometry.  The  following  boimdary 
conditions  were  applied  to  the  mesh:  (i)  for  all  nodes  on  the  line  defined  by  0=0°,  the 
displacement  in  the  direction  normal  to  this  line  was  constrained,  and  (ii)  for  aU  nodes 
on  the  line  0=2.5°,  displacement  in  the  direction  normal  to  this  line  were  constrained. 
A  typical  PAFEC  input  data  file  for  this  model  is  given  in  Appendix  A. 


uiiirillllll 


Figure  3:  Typical  finite  element  mesh  for  annulus  segment 


The  geometry  relevant  to  the  two  dimensional  idealisation  of  the  cold  expansion  of  a 
circular  hole  in  a  rectangular  plate  is  shown  in  Figure  4.  Here  the  plate  is  of  width  2w, 
length  2h,  and  has  a  centrally  located  circular  hole  of  diameter  2r.  The  analysis  of  the 
rectangular  plate  geometry  allows  a  realistic  remote  load  (i.e.  uniaxial)  to  be  applied  to 
plate  with  a  cold  expanded  hole.  Due  to  syrrunetry,  only  one  quarter  of  the  plate  was 
modelled  as  shown  in  Figures  5  and  6.  The  following  constraint  boundary  conditions 
were  appHed  to  the  mesh:  (i)  for  aU  nodes  on  the  x  axis,  displacement  in  the  y  direction 
was  constrained,  and  (ii)  for  all  nodes  on  the  y  axis,  displacement  in  the  x  direction 
was  constrained.  A  typical  PAFEC  input  data  file  for  this  model  is  given  in 
Appendix  B. 


3.2  Modelling  of  Interference  Fit  and  Cold  Expansion 

3.2.1  Elastic  Interference  Cases 

The  formulation  of  the  required  interference  constraint  conditions  in  the  finite  element 
model  can  be  described  by  considering  two  typical  adjacent  points  (A  and  B)  on  the 
mandrel/ plate  interface,  as  shown  prior  to  sleeve  insertion  in  Figure  7.  After  insertion 
these  two  points  become  coincident  for  the  no  slip  (i.e.  no  relative  displacement  in  the 
tangential  direction)  case.  Hence  the  constraint  equations  in  terms  of  Cartesian 
displacements  per  pair  of  nodes  for  this  case  are  given  in  Equations  (1)  and  (2).  These 
equations  are  then  used  for  all  pairs  of  nodes  around  the  interface.  In  the  case  of  slip 
allowed,  only  the  constraint  condition  given  in  equation  (1)  is  used.  This  enables 
relative  displacement  in  the  hoop  direction  to  occur  between  corresponding  nodes  at 
the  mandrel  and  hole  interface. 


(1) 

v„  -  -0 

p  m 

(2) 
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Figure  4:  Notation  and  geometry  for  rectangular  plate  containing  a  circular  hole 


It  shotild  be  noted  that  the  magnitude  of  the  local  cycHc  stresses  due  to  remote  loading 
are  independent  of  the  stresses  generated  by  the  applied  interference  (assuming  full 
contact  and  no  tangential  slip  at  the  mandrel/plate  interface).  Hence  the  total  stress  at 
a  given  point  in  the  plate  is  determined  simply  by  adding  the  cyclic  stress  and  the 
stress  due  to  the  interference. 

3.2.2  Plastic  Cases 

The  plastic  analyses  for  cold  expansion  and  interference  fitting  were  undertaken  in  a 
similar  manner  to  the  elastic  analyses  except  that  the  expansion  constraint  conditions 
and  remote  loading  was  applied  in  small  increments.  The  only  difficulty  that  arises  in 
using  this  method  is  specifying  correctly  the  increment  that  corresponds  to  mandrel 
removal  (ie  cold  expansion  removal)  for  geometries  other  than  the  annulus.  For  the 
plate  geometries  which  are  non-axi-symmetric  the  hole  would  end  up  non-circular 
after  cold  expansion.  During  the  iterative  finite  element  analysis,  the  increment  in  die 
cold  expansion  cycle  considered  "mandrel  removal"  was  taken  to  be  when  the  radial 
stress  at  the  critical  location  (0=0°)  becomes  zero.  The  relatively  small  error  arising 
from  this  assumption,  as  a  function  of  angular  position,  is  shown  in  Figure  8  for  the 
typical  case  of  a  cold  expanded  hole  as  considered  in  Section  5.4.  Here  the  maximum 
error  is  32.7  MPa,  and  occurs  at  0=90°.  For  comparison  purposes,  the  corresponding 
residual  hoop  and  maximum  radial  stresses  are  also  shown  in  this  figure. 
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Figure  5:  Typical  overall  finite  element  mesh  for  modelling  a  rectangular  plate  containing  a 
circular  hole 


3.2.3  General  Comments 

The  analysis  procedure  presented  above  is  considered  valid  for  the  assumed  constraint 
conditions,  with  the  further  requirement  that  no  separation  at  the  mandrel/plate 
interface  occurs  during  remote  loading  of  the  plate.  In  practice  these  requirements  will 
depend  on  interference  level,  the  relative  material  stiffness  values  of  the  mandrel  and 
plate,  the  coefficient  of  friction  at  the  mandrel/ hole  interface,  and  the  magnitude  of  the 
remote  loading. 


Figure  6:  Typical  finite  element  mesh  near  region  of  hole  boundary  for  rectangular  plate  case 


X 


Figure  7;  Generic  geometry  for  definition  of  interference  fit  boundary  conditions  along 
mandrel/plate  interface 
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Figure  8:  Stress  variation  around  hole  boundary:  for  case  presented  in  Section  5.4:  (i)  radial 
stress  at  the  maximum  cold  expansion  level,  and  (ii)  residual  radial  and  hoop 
stresses  at  completion  of  cold  expansion. 


3.3  Elastic-Plastic  Material  Properties 

Analyses  undertaken  for  D6ac  steel  (F-lllC  wing  pivot  fitting  material)  used  one  of 
the  following  two  non-linear  material  constitutive  response  curves:  (i)  an  elastic- 
perfectly  plastic  response  to  enable  direct  comparison  of  finite  element  results  with 
closed-form  theoretical  solutions,  or  (ii)  a  strain  hardening  response  to  allow  more 
realistic  simulation  of  true  material  characteristics.  These  two  curves  are  shown  in 
Figure  9.  Some  preliminary  analyses  were  also  performed  with  aluminium  plate 
material,  assuming  an  elastic-perfectly  plastic  response,  with  an  elastic  modtilus  of  69 
GPa  and  a  yield  point  of  480  MPa. 
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Figure  9:  Assumed  elastic-perfectly  plastic  and  strain  hardening  material  constitutive 
responses  ofD6ac  steel. 


4.  Benchmark  Finite  Element  Analyses 

The  following  three  plane-strain  analyses  were  undertaken  to  allow  comparison  of 

typical  finite  element  results  with  theoretical  solutions. 

1.  Elastic  analysis  of  a  large  rectangular  aluminium  plate  with  an  interference  fitting 
of  1%  due  to  a  steel  mandrel. 

2.  Elastic-perfectly  plastic  analysis  of  an  aluminium  annulus  with  an  interference 
fitting  of  4%  due  to  a  steel  mandrel. 

3.  Elastic-perfectly  plastic  analysis  of  a  steel  annulus  with  cold  expansion  to  1.5  %, 
followed  by  0.5%  interference  fitting  due  to  a  steel  mandrel. 


4.1  Elastic  Analysis  of  a  Large  Rectangular  Aluminium  Plate  with  a 
1%  Interference  Fitted  Steel  Mandrel 

For  this  analysis  the  plate  geometry  as  shown  in  Figure  4  was  considered,  with  a  hole 
diameter  of  20  mm,  plate  length  of  400  mm,  and  plate  width  of  200  mm.  The  plate 
material  was  an  aluminium  alloy  with  properties  Ep=72.4  MPa,  and  Vp=0.33.  The 
mandrel  material  was  steel  with  properties  Em=209  GPa  and  Vm=  0.33.  The  plate  was 
loaded  with  a  remote  stress  of  100  MPa  and  fitted  with  a  1%  interference  steel  mandrel. 
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The  results  for  peak  radial  and  hoop  stresses  in  the  plate  and  the  hole  edge  were 
within  2%  of  those  obtained  in  prior  work  [6]. 

4.2  Elastic-Perf  ectly  Plastic  Analysis  of  an  Aluminium  Annulus  with  a 
4%  Interference  Fitted  Steel  Mandrel 

This  case  was  chosen  to  allow  for  comparison  of  the  finite  element  results  and  the 
closed  form  solutions  presented  in  Reference  4.  The  plate  geometry  as  shown  in  Figure 
2  was  corisidered,  with  material  properties,  Ep  =  69  GPa,  Up=0.33,  ao  =  480  MPa,  and 
geometric  parameters,  a  =  10  mm,  and  b  =  100  mm,.  The  mandrel  elastic  material 
properties  were  Em  =  209  GPa  and  Om=0.3.  The  2.5°  annulus  segment  finite  element 
model  as  shown  in  Figure  3  was  used.  The  mandrel  was  represented  by  120  divisions, 
and  the  annulus  by  400  divisions.  To  model  the  interference-fitting,  the  constraint 
conditions  were  appHed  incrementally  as  follows:  one  increment  of  0.3808% 
interference  followed  by  eight  increments  of  0.0544%  interference  and  then  fourteen 
increments  of  0.136%  interference.  The  large  first  increment  was  used  to  get  up  to  just 
below  the  yield  condition,  and  then  followed  by  smaller  increments  to  get  past  the 
yield  point.  Good  correlation  between  the  finite  element  results  and  the  theoretical 
predictions  [4]  was  obtained,  as  indicated  in  Figure  10. 

The  location  of  the  elastic-plastic  boundary  at  the  4%  interference  level  is  illustrated  in 
Figure  11,  which  shows  the  von  Mises  stress  distribution  in  the  annulus  from  the  hole 
edge  to  the  annulus  outside  diameter. 


Figure  10:  Comparison  of  results  for  hoop  and  radial  stresses  for  an  interference  fitted  steel 
mandrel  in  an  aluminium  annulus  at  a  4%  interference  level 
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Figure  11:  Von-Mises  stress  distribution  as  a  function  of  radius  showing  location  of 
elastic/plastic  boundary  determined  from  finite  element  analysis  at  4% 
interference 


The  analysis  of  Jost  [5]  can  be  used  to  calculate  the  location  of  the  elastic-plastic 
boundary  (c).  For  the  given  interference  level  of  4%,  the  value  of  c  is  determined  by 
iteration  from  Equation  (3)  where 


Use  of  this  equation  yields  a  prediction  for  c  of  26.36  mm,  which  compares  favourably 
with  the  value  obtained  using  finite  element  analysis  of  24.4  mm  as  shown  in  Figure 
11.  Since  the  closed  form  solutions  [5]  are  based  on  a  deformation  theory  of  plasticity, 
and  the  finite  element  analyses  use  incremental  plasticity,  some  discrepancy  between 
the  results  is  to  be  expected. 
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4.3  Elastic-Perfectly  Plastic  Analysis  of  a  1.5%  Cold  Expanded  Steel 
Annulus  With  a  Subsequent  0.5  %  Steel  Mandrel  Interference 
Fitting 

In  this  case  the  finite  element  mesh  (Figure  3)  consisted  of  120  equal  divisions  to 
represent  the  mandrel  and  400  to  represent  the  annulus,  as  for  the  case  discussed  in 
Section  4.2.  The  geometric  and  material  parameters  used  were  as  follows;  b  =  75  mm,  a 
=  10.0  mm,  Vp  =  v„  =  0.3,  g  =  1,  oo  =  1200  MPa  and  Ep  =  Em  =  209  GPa.  To  model  the  cold 
expansion,  the  constraint  conditions  were  applied  incrementally  as  follows:  one 
increment  of  0.54%,  followed  by  two  increments  of  0.015%,  six  increments  of  0.03% 
and  ten  increment  of  0.075%.  This  resulted  in  a  total  cold  expansion  of  1.5%.  Unloading 
from  this  level  was  with  two  increment  at  -0.075%,  six  increments  at  -0.15%,  one 
increment  of  -0.075%,  and  then  one  increment  of  -0.015%.  It  is  important  to  note  that 
complete  unloading  was  defined  to  occur  when  the  radial  stress  at  the  edge  of  the  hole 
was  zero.  The  subsequent  interference  fitting  level  of  0.5%  was  achieved  using  six 
increments  of  0.075%,  one  at  0.045%  and  one  at  0.0045%. 

4.3.1  Stresses  at  1.5%  Cold  Expansion 

The  radial  and  hoop  stress  distributions  in  the  annulus  at  tire  1.5%  cold  expansion 
level,  (i.e.  before  unloading)  obtained  from  the  firute  element  analysis,  are  shown  in 
Figure  12. 


Figure  12:  Finite  element  results  for  hoop  and  radial  stress  distributions  in  a  steel  annulus  at 
1.5%  cold  expansion 
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Predicted  stresses  in  the  plastic  region  [5]  are  given  by  equation  (4),  A  comparison 
between  the  theoretical  predictions  and  finite  element  results  at  the  hole  edge  are 
given  in  Table  1,  showing  good  agreement. 


r  (  c 

21n-  +  l+  7 
c  \b 


(4) 


Table  1:  Finite  element  and  theoretical  stresses  at  the  hole  edge  for  a  steel  annulus  at  1.5% 
cold  expansion 


Component 

Theory  [5] 

Finite  element 

(MPa) 

(MPa) 

-1337 

-1287 

<^e 

48 

94 

Figure  13  shows  the  finite  element  results  for  the  von  Mises  stress  distribution  as  a 
function  of  radial  position  through  the  annulus.  This  indicates  an  elastic-plastic 
boundary  at  c  =  15.69  mm,  as  compared  to  a  value  of  c  =  16.30  mm  predicted  from 
equation  (3). 


Figure  13:  Finite  element  results  for  the  von-Mises  stress  distribution  through  a  steel 
annulus  showing  location  of  elastic-plastic  boundary  at  1.5%)  cold  expansion. 


14 


DSTO-TR-0523 


4.3.2  Stresses  After  Removal  of  1.5%  Cold  Expansion 

Figure  14  shows  the  finite  element  results  for  the  radial,  hoop  and  von  Mises  stress 
distributions  in  the  annxilus  at  the  completion  of  the  cold  expansion  process.  The 
magnitude  of  the  residual  compressive  hoop  stress  is  a  maximum  at  the  hole  edge,  and 
the  material  there  is  just  below  the  incipient  re-yield  point. 


Figure  14:  Finite  element  stress  results  for  a  steel  annulus  after  1.5%  cold  expansion  then 
removal 

From  Reference  5  theoretical  stresses  in  the  plastic  region  at  this  condition  (with  no 
reyielding  upon  mandrel  removal)  are  given  by  equation  (5).  A  comparison  between 
theoretical  predictions  and  finite  element  results  at  the  hole  edge  is  given  in  Table  2, 
showing  good  agreement.  It  is  important  to  note  that  the  residual  radial  deformation 
after  the  removal  of  cold  expansion  is  0.0366  mm,  or  0.366%  of  the  initial  hole  radius. 
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where  a  =  2  In 
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Table  2;  Finite  element  and  theoretical  stresses  at  the  hole  edge  after  unloading  from  a  1.5% 
cold  expansion  level 


Component 

Theory  [4] 

Finite  element 

(MPa) 

(MPa) 

0 

-6 

(^9 

-1338 

-1233 

4.3.3  Stress  History  for  1.5%  Cold  Expansion  and  Subsequent  0.5%  Interference 
Fitting 

Figure  15  shows  the  variation  of  hoop,  radial  and  von-mises  stresses  at  the  edge  of  the 
hole,  as  a  function  of  the  stage  (i.e.  time)  in  the  cold  expansion  process  and  the 
subsequent  interference  fitting.  The  units  on  the  horizontal  axis  represents  a  fictitious 
time  scale,  where  each  time  interval  corresponds  to  an  increment  in  the  finite  element 
analysis.  Since  the  analysis  increments  were  non-uniform,  it  is  a  non-Hnear  scale.  The 
elastic-perfectly  plastic  response  of  the  material  can  be  clearly  seen  by  inspection  of  the 
von  Mises  stress  response.  During  the  unloading  phase  of  the  cold  expansion  process, 
the  material  response  is  elastic,  with  the  von  Mises  stress  not  being  sufficient  in 
magnitude  to  achieve  re-yielding.  It  can  be  seen  that  the  insertion  of  the  interference  fit 
mandrel  induces  only  elastic  deformations  in  the  plate. 
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Figure  15 :  Finite  element  stress  history  at  the  hole  edge  of  the  steel  annulus,  for  cold 
expansion  to  1.5%,  followed  by  0.5%  interference  fitting. 


The  closed-form  solution  to  this  problem  [4]  is  given  in  equation  (6). 


CE  ^CE 

a,  =  cr,  — 


D 


-CE 


=  or, 


CE 


(6) 


Here  the  first  term  in  each  equation  is  due  to  cold  expansion  (equation  (5)).  The  second 
term  is  due  to  interference  only  and  is  assumed  to  vary  linearly  with  increasing 
interference  level.  The  results  from  the  finite  element  analysis  and  those  predicted 
from  theory  [4]  are  presented  in  Table  3,  showing  good  agreement. 
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Table  3:  Finite  element  and  theoretical  stresses  at  the  hole  edge  for  an  annulus  after  1.5%  cold 
expansion  and  subsequent  interference  fitting  to  0.5% 


Component  Theory[4]  Finite  element 


-567  MPa  -568  MPa 

-750  MPa  -652  MPa 


5.  Analysis  of  Enhanced  Circular  Hole  in  a  Large 
Rectangular  Plate  Under  Representative  F-lllC 
Loading  Conditions 

Following  the  successful  analyses  presented  in  Section  4  on  various  benchmark 
problems,  more  complex  cases  representative  of  the  F-lllC  wing  pivot  fitting  material 
constitutive  behaviour  and  loading  conditions  were  -undertaken.  For  example,  results 
of  analyses  imdertaken  assuming  a  strain  hardening  response  of  D6ac  steel  are 
presented  (Figure  9).  Also,  loads  representative  of  those  experienced  by  the  stiffener 
containing  Fl^FilS,  are  applied,  namely:  (i)  Cold  Proof  Load  Test  (CPLT)  and  (ii)  a 
sample  remote  spectrum  loading  of  +  200  MPa. 

The  plastic  finite  element  analyses  described  in  this  Section  are  for  a  large  rectangular 
flat  steel  plate  which  is  300  mm  long,  150  mm  wide  and  with  a  10  mm  radius  circular 
hole.  For  all  cases  the  finite  element  mesh  was  similar  to  that  shown  in  Figure  5,  and 
unless  otherwise  noted  the  cold  expansion  level  was  1.5%  and  the  subsequent 
interference  fit  level  was  0.5%,  and  the  no-slip  constraint  conditions  were  used 
(Equations  1  and  2).  The  results  are  presented  at  the  critical  location,  i.e.  at  0=0°  (Figure 
4),  for  both  slip  and  no-slip  cases.  The  cases  considered  can  be  summarised  as  follows: 

Unenhanced 

1.  Remote  loading  only  for  plate  with  an  open  circular  hole  assuming  strain 
hardening  material  properties. 

Enhanced 

1.  Cold  expansion  and  subsequent  interference  fitting  assuming  elastic-perfectly 
plastic  material  properties  for  the  plate. 

2.  Cold  expansion,  subsequent  interference  fitting  and  remote  loading  assuming 
elastic-perfectly  plastic  material  properties  for  the  plate. 

3.  Cold  expansion,  subsequent  interference  fitting  and  remote  loading  asstuning 
strain  hardening  material  properties  for  the  plate. 

4.  Cold  expansion,  subsequent  interference  fitting  and  remote  loading  assuming 
strain  hardening  material  properties  for  the  plate,  with  slip  allowed  to  occur 
during  the  interference  fitting  process  and  under  the  remote  loading. 

5.  Interference  fitting  and  remote  loading  for  the  plate  with  strain  hardening 
material  properties,  for  both  no-slip  and  shp  allowed  cases. 
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5.1  Unenhanced  Hole  With  Remote  Loading  Assuming  Strain 
Hardening  Material 

For  this  analysis  strain  hardening  material  properties  were  used  to  permit  direct 
comparison  with  the  results  of  the  enhanced  analyses  presented  in  Sections  5.4,5.5  & 
5.6.  The  stress  results  for  the  edge  of  the  open  hole  at  0=0°  are  given  in  Figure  16. 


Figure  16:  Finite  element  stress  history  at  the  hole  edge,  6=0°,  for  the  unenhanced 
rectangular  plate  of  strain  hardening  material  with  CPLT  and  then  remote 
spectrum  loading 


As  can  be  seen  from  Figure  16,  residual  tensile  stresses  due  to  the  CPLT  equivalent 
loading  occur,  and  this  is  consistent  with  experimental  results  obtained  from  strain 
surveys  conducted  during  F-lllC  wing  loading  tests.  The  stress  response  is  elastic  for 
the  application  and  removal  of  the  remote  tensile  loading  of  228  MPa.  However  the 
high  remote  compressive  loading  to  -556  MPa  causes  material  yielding,  and  upon  its 
removal,  a  residual  tensile  hoop  stress  of  200  MPa  is  induced.  The  stress  response  due 
to  the  subsequent  sample  remote  loading  is  Linear,  and  hence  a  residual  hoop  stress  of 
200  MPa  exists  at  the  end  of  this  loading.  Since  the  hole  is  unfilled,  the  cyclic  stress 
concentration  factor  due  to  remote  loading  is  3.1. 
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5.2  Cold  Expanded  and  Interference  Fitted  Hole  Assuming  Elastic- 
Perf ectly  Plastic  Material 

Here  cold  expansion  was  modelled  with  one  increment  of  0.21%  interference  followed 
by  eight  increments  of  0.03%  and  fourteen  increments  of  0.075%.  Cold  expansion 
unloading  (removal)  was  achieved  with  fourteen  increments  of  -0.075%  and  one 
increment  of  0.09%  to  a  achieve  a  condition  of  zero  radial  stress  at  the  critical  location. 
Interference  fitting  of  the  mandrel  was  achieved  with  10  increments  of  0.05%. 


Figure  17:  Finite  element  stress  history  at  the  hole  edge,  0=0°  for  the  enhanced  rectangular 
plate  of  elastic-perfectly  plastic  material  with  cold  expanded  to  1.5%,  followed  by 
0.5%  interference  fitting 


The  results  presented  in  Figure  17  for  this  case  serve  as  a  useful  guide  to 
tmderstanding  clearly  the  mechanism  of  stress  response  due  to  the  cold  expansion  and 
interference  processes.  During  cold  expansion,  the  compressive  radial  stress  increases 
in  magiutude  with  increasing  interference  level.  As  the  steel  was  modelled  here  as  an 
elastic-perfectly  plastic  material,  the  von  Mises  response  curve  is  flat  after  the  yield 
point  is  reached.  Hence  once  yielding  has  occurred,  the  hoop  stress  must  reduce  with 
the  increase  in  interference  (i.e.  increasing  magnitude  of  compressive  radial  stress).  On 
removal  of  the  cold  expansion,  such  that  the  radial  stress  is  zero,  there  is  a  resulting 
residual  hoop  stress  of  -1227  MPa.  The  subsequent  application  of  the  interference 
mandrel  then  results  in  residual  compressive  radial  stress,  and  a  reduction  in  the 
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magnitude  of  the  residual  compressive  hoop  stress.  The  finite  element  results  from  this 
analysis  of  the  large  plate  can  be  compared  directly  with  the  those  obtained  for  the 
annulus  shown  in  Figure  15,  along  with  the  corresponding  analytical  predictions.  For 
comparison  purposes,  the  results  at  different  stages  in  the  enhancement  process 
history  are  siunmarised  in  Tables  4,  5  &  6,  for  the  location  6=0°  at  the  hole  edge.  The 
residual  expansion  of  the  hole  boundary,  after  completion  of  the  cold  expansion 
process,  was  found  to  be  0.035740  mm  ±  0.000121  mm. 


Table  4:  Stresses  at  the  hole  edge  at  1.5%  cold  expansion  level 


Component 

Annulus 
(Theory  [1]) 

Annulus 
(Finite  element) 

Plate 

(Finite  element) 

-1337  MPa 

48  MPa 

-1287  MPa 

94  MPa 

-1212  MPa 

172  MPa 

Table  5:  Stresses  at  the  hole  edge  after  removal  of  1.5%  cold  expansion  level 

Component 

Annulus 

(Theory[l]) 

Annulus 
(Finite  element) 

Plate 

(Finite  element) 

0 

-1338  MPa 

-6  MPa 
-1233  MPa 

0 

-1227  MPa 

Table  6:  Stresses  at  the  hole  edge  after  1.5%  cold  expansion  and  subsequent  interference  fitting 
to  0.5% 

Component 

Annulus 
(Theory  [1]) 

Annulus 
(Finite  Element) 

Plate 

(Finite  Element) 

-567  MPa 
-750  MPa 

-568  MPa 
-652  MPa 

-572  MPa 
-647  MPa 

The  differences  between  analytical  and  finite  element  restdts  for  the  annulus  were 
discussed  in  Section  4.2.  From  the  results  given  in  Tables  4  to  6,  it  can  be  concluded 
that  the  finite  element  modelling  methods  used  for  the  analysis  of  the  large 
rectangular  large  plate  are  sound. 
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5.3  Cold  Expanded  and  Interference  Fitted  Hole  with  Remote  Loading 
Assuming  Elastic-Perfectly  Plastic  Material 

This  analysis  case  is  the  same  as  that  presented  in  Section  5.2,  except  that  uniaxial 
remote  loading  is  now  applied  to  the  plate  after  the  completion  of  the  interference 
fitting.  The  remote  loading  was  chosen  to  represent  the  CPLT  followed  by  a  sample 
spectrum  loading.  One  objective  for  this  analysis  is  to  provide  results  that  may  be 
compared  with  those  obtained  for  the  case  assuming  strain  hardening  material 
response,  which  are  given  in  the  next  section.  The  results  for  the  analysis  are  shown  in 
Figure  18. 


sample  remote 

0.5%  interference  fit  plug  load  spectrum 


Figure  18:  Finite  element  stress  history  at  the  hole  edge,  0=0°  for  the  enhanced  rectangular 
plate  of  elastic-perfectly  plastic  material,  with  cold  expanded  to  1.5%,  followed  by 
0.5%  interference,  CPLT  and  then  remote  spectrum  loading  (no  slip) 


Here  it  can  be  seen  that  during  the  sample  remote  loading  of  ±200  MPa,  the  mean 
stress  is  compressive  at  -647  MPa,  and  the  cyclic  stress  peaks  are  -447  MPa  and  -847 
MPa.  This  gives  a  cyclic  stress  concentration  factor  of  unity.  It  can  be  seen  that  the 
fatigue-damaging  effect  of  a  high  mean  stress  coupled  with  a  large  cyclic  stress 
concentration,  as  previously  existed  for  the  open  hole  (Figure  16)  case,  has  been  greatly 
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reduced.  Another  point  of  interest  demonstrated  by  this  analysis  is  the  constant  radial 
stress  value  during  the  application  of  remote  loading.  This  is  because  the  interference 
fit  mandrel  installation  has  effectively  filled  the  hole,  so  then  the  remote  uniaxial 
loading  does  not  create  any  stress  in  the  direction  perpendicular  to  its  application. 


5.4  Cold  Expanded  and  Interference  Fitted  Hole  with  Remote  Loading 
Assuming  Strain  Hardening  Material 

The  use  of  strain  hardening  material  properties  gives  a  more  realistic  estimate  of 
stresses  during  the  cold  expansion  and  interference  fit  process.  The  magnitude  of  the 
residual  stress  at  the  edge  of  the  hole,  and  the  size  of  the  compressive  region,  are 
significantly  influenced  by  the  work  hardening  behaviour  of  the  material.  Typically 
the  use  of  strain  hardening  material  causes  the  size  of  the  compressive  residual  stress 
region  to  be  larger  for  a  given  level  of  cold-expairsion,  as  compared  to  the  elastic- 
perfectly  plastic  case.  The  results  obtained  are  given  in  Figure  19,  and  can  be  compared 
directly  with  those  obtained  for  an  elastic-perfectly  plastic  material  as  given  in  Figure 
18. 


0.5%  interference  fit  plug 


Figure  19:  Finite  element  stress  history  at  the  hole  edge,  &=0°  for  the  enhanced  rectangular 
plate  of  strain  hardening  material,  with  cold  expanded  to  1.5%,  followed  hy  0.5% 
interference,  CPLT  and  then  remote  spectrum  loading  (no  slip) 
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It  is  seen  from  Figure  19  that  due  to  the  strain  hardening  material  response,  the 
compressive  residual  hoop  stress  (after  the  CPLT)  is  much  lower  in  magnitude,  as 
compared  to  die  elastic-perfectly  plastic  material  case,  having  changed  from  -647  MPa 
to  -1080  MPa.  The  radial  stress  remains  essentially  the  same  as  for  the  elastic-perfectly 
plastic  material  case  because  the  material  stress  response  during  the  interference  fit 
process  is  elastic  (the  radial  stress  having  started  from  zero).  The  Unear  response 
subsequent  to  interference  fitting  indicates  that  the  cycUc  stress  concentration  factor 
remains  the  same  as  for  the  elastic-perfectly  plastic  material  case,  i.e.  unity.  The 
residual  deformation  around  the  hole  after  the  cold  expansion  process  was  found  to  be 
0.030122  mm  +  0.000075  mm,  which  is  sUghtly  less  than  the  value  of  0.035740  mm 
obtained  for  the  plate  with  elastic-perfectly  plastic  material  in  Section  5.2. 


5.5  Cold  Expanded  and  Interference  Fitted  Hole  with  Remote  Loading 
Assuming  Strain  Hardening  Material  and  Slip  Allowed 

In  practice,  the  interference  fit  mandrel  can  potentially  sHp  during  the  appUcation  of 
external  loading,  if  the  loading  is  high  enough.  To  investigate  the  effect  of  sUp  on  plate 
stresses,  analyses  were  performed  for  the  case  where  there  is  no  displacement 
constraint  at  die  mandrel/ plate  interface  in  the  local  hoop  direction  (see  Equation  (2) 
and  Figtue  7).  The  results  for  this  analysis  are  shown  in  Figxue  20. 

Subsequent  to  the  CPLT,  the  mean  hoop  stresses  for  the  sUp  allowed  case  remain 
essentially  the  same  as  for  the  no-sUp  case  as  given  in  Section  5.3.  For  example  the 
value  of  mean  hoop  stress  has  only  changed  from  -1080  MPa  for  the  no-sUp  case  to 
-1071  MPa  for  the  sHp  allowed  case.  However,  it  is  seen  that  by  allowing  sUp  to  occur, 
the  cycUc  stress  concentration  at  the  hole  edge  is  increased  as  compared  to  the  no-sUp 
case.  For  example,  the  cycUc  stress  peaks  are  -720  MPa  and  -1429  MPa,  giving  an 
increase  in  the  cycUc  stress  concentration  factor  from  1.0  to  1.79,  as  compared  to  the  no- 
sUp  case.  There  is  stUl  however,  a  major  reduction  in  the  cycUc  stress  range  as 
compared  to  those  for  the  open  hole  case  (Figiue  16),  and  significantly,  the  hoop 
stresses  never  become  tensile  during  the  remote  loading. 


24 


DSTO-TR-0523 


0.5%  interference  fit  plug 


462MPa 

mean 


-539MPa 

mean 

-1071  MPa 
mean 


Figure  20:  Finite  element  stress  history  at  the  hole  edge,  0=0  °  for  the  enhanced  rectangular 
plate  of  strain  hardening  material,  with  cold  expanded  to  1.5%,  followed  by  0.5% 
interference,  CPLT  and  then  remote  spectrum  loading  with  slip  allowed 


5.6  Interference  Fitted  Hole  With  Remote  Loading  Assuming  Strain 
Hardening  Material 

The  analysis  results  for  this  enhancement  case,  assuming  no-slip,  are  shown  in  Figure 
21.  It  can  be  seen  that  the  stress  response  is  elastic  for  mandrel  insertion  and  the 
subsequent  remote  loading.  Clearly  the  mean  stresses  are  relatively  high,  however  the 
cyclic  stress  concentration  has  been  reduced  to  unity. 

The  corresponding  results  for  die  slip  allowed  case  are  presented  in  Figure  22. 
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Figure  21:  Finite  element  stress  history  at  the  hole  edge,  0=0  °  for  the  enhanced  rectangular 
plate  of  strain  hardening  material,  0.5%  interference  fitting,  CPLT  and  then 
remote  spectrum  loading  with  no  slip 


Figure  22:  Finite  element  stress  history  at  the  hole  edge,  0=0  °,  for  the  enhanced  rectangular 
plate  of  strain  hardening  material,  0.5%  interference  fitting,  CPLT  and  then 
remote  spectrum  loading  with  slip  allowed 
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Here  we  see,  as  found  in  Section  5.5,  the  main  effect  of  slippage  between  mandrel  and 
plate  is  to  increase  the  cyclic  stress  concentration  factor  from  unity  for  the  no-slip  case 
to  1.79  for  the  slip  allowed  case. 


5.7  Summary  of  Hoop  Stress  Responses  Due  to  Sample  Remote 
Loading 

In  the  previous  Sections,  5.1  to  5.6,  the  results  for  von-Mises,  hoop  and  radial  stresses 
have  been  presented  in  detail  for  various  enhancement  cases.  In  the  context  of  fatigue 
life  extension,  the  quantities  of  key  interest  are  the  mean  and  hoop  stress  responses 
due  to  remote  loading  at  the  critical  location,  0=0°.  To  allow  for  a  convenient 
comparison  of  the  effect  of  the  various  enhancement  cases,  these  resxilts  are  given  in 
Table  7.  It  is  immediately  apparent  that  there  is  a  very  significant  advantage  in  having 
combined  cold-expansion  with  interference  fitting,  as  opposed  to  interference  fitting 
only. 


Table  7:  Mean  and  cyclic  hoop  stresses  at  the  hole  edge  for  various  enhancement  cases  for  a 
rectangular  plate  due  to  sample  remote  loading  of+/-  200  MPa 


Enhancement  case  and  material 
constitutive  response  assumption 

Mean  hoop 
stress 
(MPa) 

Cyclic  hoop 
stresses 
(MPa) 

Stress 

concentration 

Open  hole -(SHM)* 

200 

-415  to  817 

3.1 

No  slip  cases 

1.5  %  CE  with  0.5%  IF  -  (EPPM)** 

-647 

-847  to  -447 

1.0 

1.5  %  CE  with  0.5%  IF  -  (SHM) 

-1080 

-1280  to  -880 

1.0 

0.5%  IF  -  (SHM) 

551 

350  to  751 

1.0 

Slip  allowed  cases 

1.5  %  CE  with  0.5%  IF  -  (SHM) 

-1071 

-1429  to  -720 

1.79 

0.5%  IF  -  (SHM) 

389 

+31  to  746 

1.79 

*  strain  hardening  material 
**  elastic-perfectiy  plastic  material 
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6.  Conclusions 


This  report  presents  the  detailed  stress  responses  for  a  large  rectangular  plate  of  D6ac 
material,  subjected  to  the  following:  (a)  enhancement  by  cold  expansion  and/or 
interference  fitting,  (b)  representative  cold  proof  test  loading  (CPLT),  and  (c)  a  sample 
spectrum  loading.  The  stress  responses  were  determined  using  elastic  plastic  firiite 
element  analysis  methods  assuming  plane-strain  conditions.  For  the  case  where  there 
is  no-slip  at  .the  mandrel/hole  interface,  and  assmning  strain  hardening  material 
constitutive  response,  (unless  otherwise  noted)  the  study  indicates  the  following: 

(i)  Results  derived  from  benchmark  analyses,  assuming  elastic-perfectly  plastic 
material  response,  compare  well  with  analytical  predictions,  and  consistent  with 
the  assumptions  taken,  provide  confidence  in  the  finite  element  procedures  used. 

(ii)  For  the  tmenhanced  open  hole,  application  of  the  CPLT  results  in  a  significant 
residual  tensile  hoop  stress  of  200  MPa,  at  the  critical  location  on  the  hole  edge. 
Hence  during  application  of  the  remote  loading,  the  mean  stress  was  200  MPa  and 
the  cyclic  stress  concentration  factor  was  3.1. 

(iii)  For  all  enhancement  cases,  the  stress  response  during  the  sample  remote  loading 
subsequent  to  the  CPLT  was  linear,  with  a  cyclic  stress  concentration  factor  of 
unity. 

(iv)  For  the  combined  cold  expansion  and  interference  fitting  enhancement  case,  the 
assumption  of  strain  hardening  material  properties  as  compared  to  elastically 
perfectly  plastic  properties,  had  no  effect  on  the  cyclic  stress  concentration  factor. 
However  as  expected,  using  strain  hardening  properties  resulted  in  a  lower  mean 
stress  ( i.e.  -1080  MPa  as  compared  to  -647  MPa). 

(v)  Enhancement  through  combined  cold  expansion  and  interference  fitting  was  a 
great  deal  better  than  interference  fitting  only.  For  the  enhanced  case  with  cold 
expansion  of  1.5%  followed  by  interference  fitting  of  0.5  %,  there  was  a  residual 
(compressive)  hoop  stress  of  -1080  MPa  as  compared  to  a  residual  (tensile)  hoop 
stress  of  350  MPa  for  tiie  interference  fitting  only  case. 

(vi)  For  aU  enhancement  cases  the  assumption  of  slip  allowed,  as  compared  to  no-sHp, 
led  to  significant  increases  in  the  cyclic  stress  concentration,  i.e.  from  1  to  1.79. 
There  was  however  much  less  effect  on  the  mean  stress.  Even  for  this  highly 
conservative  assumption,  there  was  no  tensile  component  of  hoop  stress  during 
remote  loading,  for  the  combined  cold  expansion  and  interference  fitting 
enhancement. 

(vii)  Overall  enhancement  through  combined  cold  expansion  with  interference  fitting 
is  considered  highly  beneficial,  with  the  favourable  stresses  generated  not 
adversely  affected  by  subsequent  CPLT. 
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Appendix  A 

Example  PAFEC  Input  Data  File  for  1.5%  Cold  Expanded  Annulus  With 
Subsequent  0.5%  Interference  Fitting 


c 

c 

c 

c 

c 

c 

c 

c 


FILE  :  STSTSEGPS4.DAT 

TITLE  :  INTERFERENCE  FIT,  ROUND  HOLE, 
STEEL  PLATE  AND  STEEL  PLUG, 
2.5deg  SEGMENT 

Plug  Outside  Radius  =  10.15  inm 
Annulus  Outside  Rad  =  100.0  itim 


0 

C 

C 

C 

C 

C 

CONTROL 

FULL . CONTROL 

PLANE. STRAIN 

PHASE=1 

PHASE=2 

PHASE=4 

PHASE=6 

PHASE=7 

PLASTICITY 

PHASE=9 

BASE=5555555 

CONTROL . END 

NODES 

NODE . NUMBER 
1 
2 

3 

4 

5 

6 

7 

8 

9 

10 
11 


Annulus  Inside  Rad.  =  10.0  mm 
Annulus  Material  :  Steel  E  =  209E3 

Plug  Material  :  Steel  E  =  209E3 

1.5%  Cold  expansion  %  0.5%  Interference 


AXIS.NUM 

1 

1 

1 

5 

5 

1 

5 

1 

5 

6 
6 


X 

0.0 

0.65 

10.15 

0.65 

10.15 

10.0 

10.0 

100.0 

100.0 

10.15 

10.0 


fitting 


NODES 

C  Nodes  added 
NODE. NUMBER 

for  dummy  plate 

X 

Y 

1000 

160 

2 

1001 

190 

2 

1002 

160 

30 

1003 

190 

30 

AXES 

AXISN 

RELA 

TYPE 

ANGl 
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5 

1 

1 

2.5 

6 

C 

LOCAL. DIRECTIONS 

1 

1 

1.25 

NODE 

LOCAL 

PLANE 

AXIS 

4 

5 

4 

5 

10 

C 

PAFBLOCKS 

TYPE=1 

6 

4 

6 

BLOCK  GROUP  N1 

N2 

ELEMENT 

PROP 

TOPOLOGY 

C  PLATE  BLOCKS 

1  11 

C  PLUG  BLOCKS 

2 

36210 

13 

6  7  8 

9  11 

2  2  1 

C  DUMMY  BLOCK 

3 

36210 

13 

2  4  3 

5  0 

0  0  10 

3  3  1 

C 

MESH 

1 

36210 

13 

1000 

1001 

1002  1003 

REFERENCE 

SPACING. LIST 

1 

1 

2 

400 

3 

C 

PLATES. AND. SHELLS 

120 

PLATE . NUMBER 

MATERIAL . NUMBER 

THICKNESS 

12 

13 

5.0 

13 

13 

5.0 

C 

MATERIAL 

MATERIAL . NUMBER 

E 

NU 

13 

209. 0E3 

0.3 

C  . 

. . . . Steel 

C 

PLASTIC . MATERIAL 

PLASTIC . MATERIAL 

YIELD. 

CRITERION 

UNIAXIAL . PROPS 

13 

C 

UNIAXIAL. PROPS 

1 

13 

UNIAXIAL 

TYPE 

PROPERTY 

13 

1 

1200 

,  0 

C 

YIELDING . ELEMENTS 
PLAS 

GROUP 

13 

1 

C 

Plate 

13 

2 

C 

Plug 

13 

3 

C 

Dummy 

C 

C 

CONVERGENCE 

LOAD 

MAX. ITER 

QUIT 

TOLERANCE 

1 

50 

1 

1 

C 

TOLERANCES 

REFERENCE 

TOL5 

1 

0.00001 

C 

STATE . DETERMINATION 
ALGOR 

TOL 

PATH 

2 

lE-5 

1 
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C 

INCREMENTAL 
LOAD  STEP 

1  36  112222225555555555 

*  -5  -5  -10  -10  -10  -10  -10  -10  -5  -1 

*  5555553  0.3 
RESTRAINTS 

NODE. NUMBER  PLANE  AXIS  DIRECTION 

2  2  1  2 

4  2  4  2 

C 

RESTRAINTS 

NODE. NUMBER  PLANE  DIRECTION 

1000  2  0 

GENERALISED. CONSTRAINTS 
DIRECTION=l 
node  LIST 

7  5,1,1  1002,1,1 

6  3,1,1  1002,1,1 

11  10,1,1  1002,1,1 

C 

DISPLACEMENTS . PRESCRIBED 

NODE. NUMBER  DIRECTION  DISPLACEMENT . VALUE 

1002  1  0.15 

C 

END . OF . DATA 
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Appendix  B 


Example  PAFEC  Input  Data  File  for  Rectangular  plate  with  Cold 
Expanded  and  Interference  Fitted  Hole  with  Remote  Loading 

C  File  :  STSTAPS2.DAT 

C 

C  PAFEC  8.1  DATA  FILE 

C 

C  TITLE  :COLD  WORKING  1.5%  +  INTERFERENCE  FIT  0.5% 

C  CPLT  SPECTRUM  +  TYPICAL  TENSION  +  TYPICAL  COMPRESSION, 

C  ROUND  HOLE, 

C  STEEL  PLATE  AND  STEEL  SLEEVE  ;  E1/E2  =  1 

C  PLANE  STRESS  ANALYSIS 

C 

C  Quarter  plate  model  75  mm  x  150  mm 

C  Hole  radius  =  10  mm 

C  Mandrel  radius  =  10.1  mm 

C  Plate  and  mandrel  steel,  E=  209E3 

C 


CONTROL 

FULL . CONTROL 

PLANE. STRAIN 

PHASE=1 

PHASE=2 

PHASE=4 

PHASE=6 

PHASE=7 

PLASTICITY 

PHASE=9 

BASE=5555555 

STOP 

CONTROL . END 
NODES 

NODE . NUMBER 

X 

Y 

1 

0 

10.0 

2 

3.82683 

9.23880 

3 

7.07107 

7.07107 

4 

9.23880 

3.82683 

5 

10.0 

0.0 

6 

75.0 

0.0 

7 

0.0 

75.0 

8 

75.0 

75.0 

9 

0.0 

150.0 

10 

75.0 

150.0 

11 

0.0 

0.65 

12 

0.2487442 

0.6005217 

13 

0.4596194 

0.4596194 

14 

0.6005217 

0.2487442 

15 

0.65 

0.0 

16 

0.0 

10.1 

17 

3.865103 

9.331183 

18 

7.141778 

7.141778 

19 

9.331183 

3.865103 

20 

10.1 

0 

C 
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C 

C  Nodes  added 
1000 
1001 
1002 
1003 
C 

C  Nodes  added 
21 
22 

23 

24 

25 

26 

27 

28 

214 

215 

216 

217 

218 

219 

220 
221 

585 

586 

587 

588 

589 

590 

591 
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